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Abstract 


A numerical procedure is presented for computing the unsteady transonic 
flowfield about three-dimensional swept wings undergoing general time 
dependent motions. The outer inviscid portion of the flow is assumed to be 
governed by the modified unsteady transonic small disturbance potential 
equation which is integrated in the time domain by means of an efficient 
alternating direction implicit approximate factorization algorithm. Gross 
dominant effects of the shock-boundary- layer interaction are accounted for by 
a simple empirically defined model. Viscous flow regions adjacent to the wing 
surface and in the trailing wake are described by a set of integral equations 
appropriate for compressible turbulent shear layers. The two-dimensional 
boundary-layer equations are applied quasi -statically stripwise across the 
span. Coupling with the outer inviscid flow is implemented through use of the 
displacement thickness concept within the limitations of small disturbance 
theory. Validity of the assumptions underlying the method is established by 
comparison with experimental data for the flow about a high aspect ratio 
transport wing having an advanced airfoil section. 


Draft Report received August, 1982, Revised Final Report received 
November, 1983. 
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NOMENCLATURE 


a 

A 

AR 

b 


c 

c„ 


(C )rn 
t EQ, 


dTT 

m. 


V 


D 

n 


■ 'I 1 * 2 G S 5y *„ 

B e.M 2 k 2 

1 00 

* aspect ratio 

■ * S S 
3 2 M 2 k 

CO 

local section chord nondimensionalized by c r (c /c ). 

3 reference chord 
“ local section chord 

= entrainment coefficient (rV-J- / 6 P U dN) 

p e e 5 o 

= skin friction coefficient (see Eq. A-ll) 

■ total wing lift coefficient <2 j jJ , C J . Cp *>d 5 d„) 

■ total wing moment coefficient, moment taken about leading edge 
of root section 

= local section pressure coefficient 
= see Eq. A-16 

* see Eq. A-20 

* see Eq. A-13 

= type-dependent finite difference operators 
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= 1 - M c 


=> local Instantaneous surface section definition {z = f { x , t) ) 

* local viscous ramp surface slope modification 

■ - 1/2 (y + 1) M J 
- see Eq. A-15 
= see Eq. 24 
a see Eq. 25 

= Initial $ distribution, , q, 0) 

* 1/2 (y - 3) M 2 


* 1 - M 


= 1/2 (y - 1) M J -1 

= initial (j> t distribution, {£, n. 0) 
= - (y - 1) M 0 2 

§1 

= shape factor, 0 


= see Eq. A-14 


® shape parameter, 0 J Q Pg (1 - U e ) dN 
= nondimensional time scaling 


= reduced frequency, TJ 


= Mach number 

= physical normal boundary-layer coordinate 
= turbulent Prandtl number 
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adiabatic recovery factor, ( Pr^) J,/ ' 5 

p U c„ 

<a r 

freestream Reynolds number based upon chord, y 


5 c V 5 « 


P e V 

“ momentum defect thickness Reynolds number, 

3 physical streamwise boundary-layer coordinate 
3 Sutherland law gas constant 

c 

3 physical time nondimensional ized by -g- 

k ra 

3 temperature 

3 streamwise velocity component 

3 physical streamwise coordinate nondimensional ized by c f 
3 physical spanwise coordinate nondimensional ized by c 

r 

3 physical normal coordinate nondimenslcnalized by c 

r 

3 angle of attack 
3 viscous ramp wedge angle 
3 specific heat ratio 
3 jump in pptential across wake 
3 boundary-layer thickness 

pU 

3 displacement thickness, /“ (l-pIT)dN 

0 6 C 

3 finite difference operators 


A 5 »An»A(;,At = finite difference numerical mesh step sizes 


E l* e 2 


0 low frequency approximation 

1 otherwise 


= nondimensional transformed normal coordinate 
3 nondimensional transformed spanwise cordinate 


5 




\E 

A s 

X 


K 

*o 

Sp 

*s 

p 

<t> 


a gU_ U 

■ momentum defect thickness, J Q p g U e (1 - U g ) dN 

= see Eq. A-18 

3 leading edge sweep angle 
= local section sweep angle 
3 taper ratio 
3 see Eq. A- 17 
3 viscosity coefficient 

3 nondimensional transformed streamwise coordinate 
3 offset of viscous ramp leading edge from sonic point 

3 length of viscous ramp precursor 

3 length of viscous ramp main body 

3 sonic point location 
= density 

3 perturbation velocity potential function nondimensionalized by 


0 3 physical velocity potential function, c r U o (x + <*>) 

u 3 circular oscillation frequency 



3 upper surface 


a lower surface 


3 backward difference 


3 first intermediate solution (prediction) 

3 second intermediate solution (first correction) 
3 critical value 
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i 

i 

$ 


A 

■ symbolic representation 

n 

■ previous time level 

n + l 

■ final solution at new time level (second correction) 

Subscripts 


e 

* boundary-layer edge 

I. 0. K 

® n» t grid Indices 

K P 

“ 5 - plane lying above wing surface 

K m 

3 5 - plane lying below wing surface 

LE 

a leading edge 

min 

= minimum 

tip 

= wing tip 

t 

3 time derivative 

US 

■ upstream 

x,y,z 

S.n.c 

= spatial derivative 

2D 

3 two-dimensional 

3 D 

3 three-dimensional 

O 

3 freestream 
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I. INTRODUCTION 


Solutions of steady inviscid transonic flowfields about wings of arbitrary 
planform are commonly obtained by numerical solution of the differential 
equation governing the velocity potential function. While this approach 
necessarily Implies that any embedded shock waves are sufficiently weak that 
rotational ity effects may be neglected, it has proven extremely useful for 
many applications 1- ^. In the case of unsteady flows, few methods capable of 
generating time accurate three-dimensional flowfields are currently 
available . If the wing may be considered thin, the problem is rendered 
more tractable by applying the small disturbance assumption. This not only 
simplifies the governing equation, but also facilitates application of the 
surface boundary condition. 

Time integration of the modified unsteady transonic small disturbance 
potential equation has proven to be particularly useful. It permits the 
treatment of nonlinear unsteady three-dimensional flow phenomena including 
Irregular shock wave motion. In addition, the flowfield equation may be 
coupled directly with the structural equations of motion for an elastic wing 

which are simultaneously integrated in time in order to obtain a nonlinear 
13 

aeroelastic solution . 

Although unsteady inviscid flowfield solutions can provide a reasonable 
physical description for a wide variety of flow conditions, they will not be 
adequate when viscous effects are significant. This is typically true for 
flows about wings having advanced airfoil sections which involve complex flow 
structures characterized by the following phenomena: 1) shock-boundary-layer 
Interaction altering the shock strength and location, 2) effective camber 
modification due to differences in the boundary-layer displacement on the 
upper and lower wing surfaces, and 3) displacement and camber effects of the 
near wake. A detailed description of these effects can be provided only 
through the use of additional equations governing the flow. Numerical 
solution of the more exact equations appropriate to three-dimensional viscous 
flows is burdensome in terms of computing resources. For the study of 
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aeroelastlc behavior and nonlinear flutter, It Is reasonable to perform a more 
simplified analysis. It Is the intent here to account for the gross dominant 
effects of Ir.viscl d/viscous Interaction for such applications without 
degrading the basic efficiency of the inviscid computation. Thus, no attempt 
will be made to resolve exact details of viscous regions, but rather the 
effects of viscosity upon the unsteady surface pressure distribution will be 
determined In order to provide appropriate aerodynamic forces for structural 
applications. 

It is assumed that viscous effects are confined to thin regions immediately 
adjacent to the wing surface and along the trailing wake. These are then 
postulated to be in instantaneous equilibrium with the unsteady outer inviscid 
flow. A simple order of magnitude analysis indicates that this assumption is 
valid if the reduced frequency of the unsteady wing surface motion is 
small*^. As a consequence, a steady form of the equations governing the 
viscous regions may be applied in a quasi-steady fashion. In addition, the 
viscous equations appropriate for two-dimensional flows will be applied 
stripwise along the span. For this calculation, however, the streamwise 
direction is taken as that either normal to the mean of the leading and 
trailing edge sweep angles for locally subsonic strips, or normal to the shock 
surface for strips where shock waves are present. Thus, the treatment of the 
viscous equations embodies the concept of simple sweep theory. 

The use of two-dimensional boundary-layer equations greatly reduces the amount 
of computational effort involved in treating viscous regions. Even with an 
Integral formulation, a three-dimensional viscous analysis requires solution 
of a system of partial differential equations, which can be costly. This is 
no small consideration for performing aeroelastic analysis which may require a 
large number of solutions to establish flutter boundaries. 

Viscous regions are assumed to be governed by a set of boundary-layer 

equations appropriate for turbulent flows. For this purpose an integral form 

of "lag entrainment equations" are employed, which have been used with success 

15-17 

for the prediction of turbulent shear layers . Coupling with the outer 
Inviscid flow is provided by the displacement thickness concept within the 
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framework of small disturbance theory. The impact of viscous effects upon the 

potential equation is manifested solely in a modification of the surface and 

wake boundary conditions. This procedure has proven adequate for a number of 

18-20 

both steady and unsteady computations of flows about airfoils 

The most complex aspect of viscous transonic flows over airfoils is that of 
the shock-boundary-layer interaction. This effect is difficult to represent 
and resolve correctly for turbulent flows, and troublesome to couple in a 
numerically stable fashion with the outer inviscid solution. In the case of 
steady flows, procedures have been developed which achieve this coupling. 

Such methods commonly involve the use of underrelaxation in an iterative 
process such that a matching of the two solutions is obtained at convergence. 
Both the physical and numerical difficulties are considerably more involved in 
the unsteady case. Here the shock strength and location vary in time as the 
shock traverses the wing surface. In addition, time accuracy as well as 
stability is essential. 

Much of this difficulty is removed by employing a simple computational 

artifice to model the displacement effect of the shock-boundary-layer 

Interaction which is used in conjunction with the boundary- layer solution. 

This is accomplished by augmenting the wing surface geometry with a 

wedge-nosed ramp that is inserted near the base of the shock in the inviscid 

calculation. The ramp converts the normal shock to an oblique shock, thereby 

decreasing its strength, displaces it upstream from its inviscid location, and 

Is free to traverse the wing surface in correspondence with the unsteady shock 

motion. The magnitude of the ramp angle is allowed to adjust instantaneously 

to the shock strength and is chosen such that an empirically defined 

post-shock pressure is recovered. Application of this ramp model has proven 

quite useful in both two- and three-dimensional steady computations and 

14 2 2 

more recently in two-dimensional unsteady solutions ’ . 


10 


II. UNSTEADY POTENTIAL EQUATION 


The equation governing the outer inviscld flow is the modified 
three-dimensional unsteady transonic small disturbance potential equation 
which includes higher order terms to ensure proper swept shock jump 
conditions, and may be written as 

■ T % < E *x * F *x 2 * S 2) 


Here cj» is 

+ gy (<i>y + H$ x ) + g z (<{> z ) . 

the perturbation velocity potential function and 

(1) 

A 

= fc 2 , 

(2a) 

B 

- 2 M co 2 k ’ 

(2b) 

E 

CVJ 

s 8 

H 

It 

(2c) 

F 

= - 1/2 ( Y + 1) M 2 , 

(2d) 

G 

- 1/2 (y - 3) M o 2 , 

(2e) 

H 

= - ( Y - 1) M 2 , 

1 CO * 

(2f) 


where k is the nondimensional time scaling. 


The following shearing transformation is applied in order to map the swept 
tapered wing planform in the physical domain into a rectangle in the Cartesi 
computational domain: 


x - x LE (y) 

— cTyl » n = y» q = z. 


( 3 ) 
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Outboard of the wing tip, a smootn transition to a constant chord section 
swent at the average of the leading and trailing edge sweep angles Is 
applied. In addition, the coefficient G is split into two components 


G ■> G s + G n , 

(4a) 

G = 1 - M 2 , 

o o 

(4d) 

G„ = 1/2 (y - 1) -1, 

(4c) 

such that the streamwise and normal contributions to the spatial differencing 
of the governing equation may be identified. This decomposition will be 
useful in the formulation of the numerical solution algorithm. Upon defining 

a = ( !> + 2 g s 

(5a) 

b = ^2* + G S 5y 2 * 

(5b) 

X = G N { h + * n )Z + S ^ ( ?y \ + «}> n ) 

(5c) 

+ c * y Uy 


Y = ( Cy ^ + <b n ) + CE y ^ 

(5d) 

Eq. 1 becomes 


cA *tt * V k la * b ‘t 2 ’ * 2 G s » 5n 
♦§r( c * ) . +15 *f . 

(6) 

It Is noted that 


3? (a + b + 2 G s * n % 

(7) 


represents the streamwise contribution to the governing equation with the 
remaining terms on the right-hand side of Eq. 6 comprising the normal portion. 


12 



A formal definition of the Inviscid problem Is completed by prescription of 
the boundary and Initial conditions: 


far upstream $ = 0, (8a) 
far downstream ^ <j)^ + e 2 k 4> t “ 0, (8b) 
wing root 5y + = °» (8c) 
far spanwise = 0, (8d) 

far above and below $ = 0. (8e) 


On the wing surface, the linearized unsteady boundary condition 


♦ 1 + . . + 
V = cf, + c ' kf * 


: l* T t 


( 9 ) 


is appl led on <; = 0* for 0 < 5 < 1, 0 £1 1 n^p> 


Along the trailing vortex sheet in the wake, the following contact conditions 
are Invoked: 


+ e 2 k 4) t ] = °, (10a) 

UjJ 3 0, (lOo) 

on c « 0 for 5 > 1, 0 < n <. n t f p t where brackets denote the jump in the enclosed 
quantity from above to below the vortex sheet. 

Finally, the initial conditions 

0 (C.n.CiO) = 9 (?.n.c5. $ t (?.n»i;» 0 ) = <1(5, n , c ) 

(ID 


are specified. 

Once the solution to Eq. 6 is obtained, the local instantaneous pressure 
coefficient may be determined from the following expression 

C p = * 2 4*5 + e l k *t J - (12) 
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III. VISCOUS EQUATIONS 


Solution of the equations governing viscous flow regions Is predicated upon 
the quasi-study assumption. Therefore, It is necessary to consider only the 
steady form of the equations Involved In order to extend the flowfield 
solution from time level t n to time level t n+ * = t n + At during the 
unsteady solution process. In order to evaluate viscous parameters at the new 
time level, t n+ *, all Inviscld quantities are obtained from the potential 
solution at time level t n . The two-dimensional form of the viscous 
equations is employed stripwlse along the span. However, the streamwise 
direction Is taken as that normal to the local sweep angle, a $ . The value 
of A $ is equal to that of the incident shock surface for locally 
supersonic strips or the mean of the leading and trailing edge sweep angles 
for locally subsonic strips. 


Viscous Ramp 

The viscous ramp model, which is used to simulate the displacement effect of 
the shock-boundary-layer interaction, is based upon the observation in many 
steady experimental measurements that the post-shock pressure level for 
turbulent flow over an airfoil corresponds approximately to that of the 
oblique shock produced by flow over a compression ramp with a wedge angle 
equal to that of the detachment angle for the given upstream Mach number. 
Geometry of the model is depicted schematically In Figure 1. It consists of a 
short precursor over which the ramp slope varies linearly from zero to the 
given wedge angle followed by the main ramp body along which the slope varies 
quadratically. At the ramp leading edge the ramp height and slope are 
continuous, and at the downstream end and the slope and curvature vanish. The 
ramp is positioned with respect to the local instantaneous sonic point 
location and affinely scaled with the shock strength as determined by local 
conditions upstream of the shock. By offsetting the ramp leading edge a 
distance Q 0 ahead of the sonic location, the ramp is able to more fully 
influence the numerical shock profile. The leading edge of the ramp is 
preceded by a precurser of length 5 which has been used in previous 

r 
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j i a pH 22 

unsteady calculations ’ * to moderate passage of the sharp leading edge 

I across computational mesh points, thereby precluding spurious numerical 

Instabilities. The main body of the ramp has a length of £ R where the 
| parameters £ p » E 0> and may be selected for each specific 

| application. 

I • 

) . 

I 

The wedge angle, 0, Is obtained directly from the steady two-dimensional 
form of Eq. 6. A complete derivation may be found In Ref. 20. In addition, 
the concepts of simple sweep theory are applied to obtain 


6 a (P 


[ 


-2 cos 2 a. (£) <J> r - E 
s c * US 


3/2 

■] 


(13) 


where <*,_ Is evaluated just upstream of the shock. Given 4, and 

^ US * US 

the sonic point location, r , the ramp slope variation f_ Is 

K 

defined as follows: 


f * - 0 for c < 5 S - 5„ - S p . 

< * . [i . lililli"'] for 
K C ^R 

?s ' ^0 " ^pi^i^s " V 


(14a) 


(14b) 


f * ■ * U - <1 for 
\ '’R 

- 5 0 ^iq^o + ?R 


(14c) 


f R - = 0 for 5 > 5* - C 0 ♦ 5 R . 


(14d) 


IS 


£ 



Lag Entrainment Equations 


The form of the Integral boundary-layer equations which will be considered 

15 

here Is that of the lag entrainment equations due to Green . They are 

predicated upon the boundary-layer assumption that the normal extent of 

viscous regions Is small compared to the wing or wake thickness, which 

necessarily applies to flows at high Reynolds numbers. By Integrating the 

governing partial differential equations in the normal direction and suitably 

modelling the requisite relationships among the dependent variables, a system 

of three first order ordinary differential equations is obtained. Two of the 

equations result directly from continuity and momentum. The third evolves 

23 

from the Bradshaw, Ferriss, and Atwell turbulent energy equation, but is 
formulated in terms of the entrainment concept originally proposed by 
Heac£ 4 . This yields a streamwise rate equation governing the degree to 
which the outer inviscld flow merges with the turbulent shear layer. 

The lag entrainment equations are beiefly summarized here, where the simple 
sweep concept has been employed for coupling with the inviscid solution. We 
first define 


displacement thickness 6* 



(1 


P j| ) <JN, 

p e u e 


y y 

momentum defect thickness 0 - Jo p5„ (1 - S > dN ’ 

p e u e e 

jrk 

shape factor H Q = g— , 

shape parmeter IT = -£■ (1 - jp) dN, 

entrainment coefficient C r = — pU dN. 

E p e U e dS J o M 


( 15 ) 

( 16 ) 

( 17 ) 

( 18 ) 
( 19 ) 


Here S and N are the streamwise and normal directions respectively within the 
boundary layer. The subscript e refers to the boundary- layer edge denoted by 
N = o. Within the confines of small disturbance theory the edge conditions 
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may be taken as those corresponding to the tnviscfd solution along the wing 
surface or wake centerline evaluated on 5 ■ 0 *. 


0 — 

The-prlmary dependent variables are taken as ( 7 -). H, and C F . Given their 

c s t 

value at any streamrtse station, thefr distributions may be predicted by 
the following system of first order ordinary differential equations: 

„ C* cos A. 


d <c?> = ~ 2 ” " C«o + 2 - H e 2 cos2 A s ] ( c7 } hV (20) 

c f cos A s 


dtf 


► 

/ G \ — 1 | 

's' ™il 


C E cos A $ - Hj [-i— 2 - (H q + 1) ( 


(21) 

TC V ( cf rl F A { ( H7 I H7 ) cos A s I ' (C t ) EQ 2 " A c (C t ) 1/2] + ( u| "df ) EQ cos A s 


? « 1+0.2 M 2 cos 2 A. . , 

- [1 + 0.075 H 2 cos 2 A ( S- -_!)] (-i) (I) a 

e s 1 + 0.1 H / cos 2 A. c s c K 


] 


( 22 ) 


A definition of the parameters necessary to complete the description of these 
equations Is given In Appendix A. 


Using Eqs. 17 and A-12 it may be shown that 


<fj>S ■ F 1 * f 2 


(23) 


where 

F i ■ ( *-p- * 11 ♦ H e z J ic e - ^ ] cos A s (2,) 

and 

F 2 =|[1 + (3LJ_L)r M e 2 ] (H 0 + 1) Hj ^ + (y - l)r M fi 2 il + (I^I)r M g 2 ] (TT + 1) 

- H 0 (H q + 2 - He 2 ) J (f-) (i). (25) 

This form was first deduced by East^ et al. and displays the dominant 
explicit dependence of the displacement surface slope upon the pressure 
gradient <j>^. 
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Modified Boundary Conditions 


The effect of the viscous flow regions Is accounted for in the invlscid 
solution by modifying the surface boundary condition to incorporate viscous 
displacement. Thus, Eq. 9 Is replaced by 


+ 


cf„ 


+ Cjkf^* + 



where the viscous ramp is used, and 




(26) 


(27) 


on c = 0* for 0 < 5 < l, o < n < „ t1 

when coupled with the lag entrainment equations. Downstream of the wing 
trailing edge the viscous wake generates an effective displacement afterbody. 
Due to its presence, a discontinuity In the slope of the potential is now 
permitted such that the wake condition Eq. 10b is replaced by 

C V = 3 cf 1 ♦ F 2 4^] 


on 5 - 0 for 5 > l, 0 < „ < „ 


(28) 
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IV. NUMERICAL METHOD 


Invlscld Algorithm 

The time accurate solution of the Invlscld potential equation (Eq. 6 ) is 
obtained by the following first-order (in time) accurate approximate 
factorization alternating direction implicit (ADI) algorithm: 


1 ) £ - sweep: 


B 6 r ( 


<i> - <t> 


.+r + 


X ‘-At - 5 = d e Ca " * b ♦*" ^ 

+ 2 G s (6^ *") *") + 6 n (c 6 n 4 » n ) 

+ c 6 4." + 6_ X n + 6 Y n , 

•.v l n 


11 ) q - sweep: 


b 5 ? \ (c «„ * - C s n ,") 

* G s <«„ »"> » 0 <* c 4 - < r * n l. 


ill) 5 - sweep: 
cA ( 


>" + 1 - 26 ° + *" - 1 


26" + 6" ' \ » "■ ,* n + 1 . I 

t 2 2 ) + B 6 (2 1-2 


At 4 


At 


(29a) 


(29b) 


„ £ # n + 1 n. 


(29c) 


Here, and are second-order accurate central difference 

operators, ^ Is a first-order accurate backward difference operator, and 

D r and D n are type-dependent mixed difference operators. Details of 

the numerical method are summarized In Ref. 10 and discussed In detail in Ref. 
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Viscous Equations 


The viscous ramp slope modification at time level n+1 { f p + *) Is evaluated 

based upon the potential distribution at time level n ((}> n ). In order to 
locate the sonic point, the simple sweep relationship 

C_ ■ C 0 /cos 2 a (30) 

P20 P 3D s 

Is employed, where the steady form of C D Is taken as is consistent with the 

k 3D 

quasi-steady assumption. The sonic point Is then defined such that 

C = C at r = r . (31) 

P2D P20 * 

★ 

The critical pressure coefficient, C , is that corresponding to the 

P 2D 

equivalent two-dimensional freestream Mach number 

= H cos A.. (32) 

2D “3D S 

The upstream location at which Is evaluated corresponds to that grid 

‘'US 

point lying two mesh point upstream of £ . Once r and $ are 

s s HJS 

obtained, the entire ramp geometry is defined by Eqs. 13 and 14. 


Solutions to the lag entrainment equations in the form given by Eqs. 20 to 22 
are obtained by a simple second-order accurate Runge-Kutta Integration 
scheme. Initial conditions are established upstream of shock waves for 
sections where the flow is locally supersonic or at the point of minimum 
pressure for sections where shock waves are absent. For this purpose 

Incompressible flat plate results and a 1/7 power law velocity profile^ are 
used to initialize e, F, and C|r. Integration is then carried out to the 
downstream boundary based upon the potential distribution at the previous time 
Step. 
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Coupling Procedure 


Far upstream of shock waves all viscous effects upon the Invlscid solution are 
neglected. For locally supersonic stations the viscous ramp Is employed ahead 
of the sonic location as a simple wing surface slope modification given by 
Eq. 26, but no coupling with the boundary-layer solution Is involved. 
Downstream of the sonic location the remaining portion of the viscous ramp Is 
ignored and the boundary-layer solution is implicitly coupled to the Inviscid 
solution through use of Eqs. 27 and 28. In the case of locally subsonic 
stations no viscous ramp is employed and the Implicit coupling procedure is 
Initiated at the point of minimum pressure. 

The wing plane = 0) is centered between planes of computational mesh 
points, such that the surface boundary condition is manifested In the 
calculation of ^ In the mesh planes just above and below ^ = 0. This 
is accomplished by employing the relationship 


^IJK +1 ‘ *IJK 5 " ^ (33) 

P P P 

where the index K p refers to the mesh plane just above the wing surface and 
a constant mesh size is considered here for illustrative purposes. The 
surface boundary condition is now employed as 

\ } ij = cj ^IJ + c l k(f t l IJ + (F I + F 2 * 34) 

and is written in terms of the potential at points lying 

above the surface; i.e. 


(<j) E5 ’iJ = ( a| ) Cc K +1 \ +1 ] * 

r r K r 
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Equations 34 and 35 are then used In conjunction with Eq. 33 to obtain 

/ %n + 1 # lv2 f , n + 1 ,n + 1. 

'♦cc'IJK 'Ic' | ^IJK + 1 ^IJK ' 

r r r 

. r 1 jr +*n + 1 +%n + 1 . /c +%n n 

■ ^ ’■Cj (f E *IJ e i k(f t 1 IJ (F l *10 3 

_ ( + i n r_ (x^ _o x** + I + x n + ^ ) 

1 2 'IJ L? K + l l + I + 1 JK 1 't’lJK *1 - 1 JK 
s P p p p 

“ c K p (<!, i + i jKp+ i " 2 *uy 1 + ♦I - 1 J K p + l*- 1 ] • (36) 

This expression forms an implicit relationship as part of the tri-diagonal 

system in the x, - sweep of the invlscid algorithm. An expression similar to 

Eq. 36 may be derived for (<!>,., )ijk *• It is noted that at the grid location I, J valu 

of 4 ) at the upstream column (♦” t J j) are known at the advanced time level but 
those at the downstream column (4>i + l j) must be taken at the previous time 

level for the algorithm to remain intact. The coupling procedure is considered 

to be significant in suppressing numerical Instabilities which would otherwise 

★ 

X 

occur, particularly in regions where (£— ) r is large such at the trailing edge. 

c s * 

Along the trailing wake the following relationship is employed 


^ 55 * IJK J (<!> IJK +1 ’ 2<!> IJK + *IJK 1 (37) 

p p pm 

where is a symbolic representation of the value of a t ,„ which 

idK m IdK m 

incorporates the jump in $ and <|>^ across 5=0 and K m refers to the mesh 
plane just below 5 = 0 . By analytic continuation (i.e., Taylor Series) 
4 >j JK ^ may be expressed in terms of values along 5 = Q + , i.e. 


= ‘hj “ ( r _) ( Vu. 
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Similarly, the value of <{>jj|, may be written In terms of values along ^ = 0" as 


*lJK m = *IJ ‘ {< t>^U • (39) 

Combining Eqs. 38 and 39 results In 

*IJK m ’ ♦uk, * <*u - <*> - (^1 t(**)„ - (♦')„]. (4()) 

For convenience we now define 

{ *IJ ' ♦iJ* = r u (41) 

and express the jump in the normal slope of the potential across the trailing 
wake in accordance with Eq. 28 as 

C( *C , IJ ' ( *C } I J ] = ilF l ] U + (F 2 } IJ 

" ^ F 1 J IJ + (F 2 ) IJ (42) 

From Eq. 41 It follows that 

(*«)lJ = (43) 

Using Eq. 35 and substituting Eqs. 41 to 43 into Eq. 40, it may be shown that 

*UK m ’ *IJK m + r XJ - (£) C(Fl) u - (F p H ] 


- 7 - ^> U 3 tty! V, JKp • \ +1 ] 


" ‘j 5 - 1 (F 2>M (r rr l lg- 
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Finally, Eq. 44 Is substituted Into Eq. 37 to obtain 

i , »n+l _ i 1 \2 # n+1 «.n+l . .n+1 + ~n+l 

( *CC UKp " W ^IJKp+1 Z ^I0K p + r IJ 


" [(F 1 1 IJ ’ (F 1 , I0 3 


(i) a) 2 a> 2 C(F*)" U - (r 2 )- u ] eck +1 <»W - < ♦ OK > 


AC* 4E 


P 

" ^K p ^I+l JKp+1 _ 2 $IJK p +l + *1-1 JKp+1^ 


- <7S;> ( i |2 <9u <U 


2 n+1 n+1 . 

zr I0 + r l-lJ J * 


(45) 


n+1 


An expession similar to Eq. 45 may be derived for Uj-^ijk 


f 

l 

? 


! 

i 
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V. RESULTS 


A computer program, XTRAM3S, based upon the Inviscid numerical algorithm 
described above is available for calculating the steady and unsteady 
aerodynamic loads and aeroelastic response of thin clean wings in transonic 
flow. This code is fully described in Ref. 25 to 28. All of the results 
presented here were generated using XTRAM3S. The basic code has been modified 
to account for viscous effects by incorporating the previously described 
method. Appendix B describes alterations in user input specifications for the 
modified version of XTRAN3S. 

For each of the cases considered, calculations were performed which generated 
an inviscid colution, a solution using the viscous ramp alone, and a full 
viscous solution employing the viscous ramp in conjunction with the 
boundary-layer equations. 

All computations were performed on a nonuniform 60 x 20 x 40 (£, n, 

Cartesian computational mesh with the wing surface defined by 39 x 12 points 
In the E - n plane. The computational domain was defined by 

-15.375 < C <26.575 
0 < n 1 5.3 

-13.0375 < c < 13.0375 


and minimum grid spacing taken as 



which occur at the wing leading edge, at the wing tip, and adjacent to the 
wing surface respectively. XTRAN3S default values were employed for the 5- 
and c-mesh distributions. The n-mesh and surface geometry description for 
each case considered may be found in Appendix C. For all calculations the 
time scaling, k, was selected as 0.2 and a time step of At = 0.0034906585 
was employed. This corresponds to a distance of one root chord of travel in 
57.3 time steps at the freestream velocity. For a reduced frequency of 
k £ = 0.2, the choice of At results in five time steps per degree of 
circular frequency change along the pitching cycle of a forced oscillation. 
Nominal values of the viscous ramp parameters were selected as follows: 

e 0 « 0.02, 

K p = 0 . 02 , 

5 r 3 0 . 10 . 

These choices have proven adequate for a number of both steady and unsteady 
two-dimensional solutions 14,20 . Results were generated on the CRAY-1S 
computing system and required approximately 0.5 seconds of CPU time per time 
step of calculation for both inviscid and weage alone solutions, and Q.b 
seconds per time step for full viscous computations. 

All converged steady state solutions were run for 900 time steps at the 
indicated values of k and At. These choices were found to be conservative 
with respect to both stability and convergence. The low frequency 
approximation (i.e., e 2 3 °) was made <n the wake Jump condition Eq. 10a 
and downstream boundary condition Eq. 8b to prevent anomalous behavior which 
was otherwise observed to occur for swept tapered planforms. In general, the 
following procedure was employed in generating steady-state results: 

1) a converged inviscid solution was obtained using an undisturbed 
condition as an initial state {i.e., $ * $ = 0); 

t 

ii) a converged wedge alone solution was obtained using the inviscid 
solution as the intiial state, updating the wedge computation at 
each time step; 
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Ill a converged full viscous solution was obtained using the wedge 
alone solution as the initial state, updating the boundary-layer 
computation at each time step. 

In order to establish validity of the assumptions underlying the computational 

method, comparison is made with experimental data for the steady flow about a 

high aspect ratio wing having an advanced airfoil section. The configuration 

to be considered Is depicted schematically in Fig. 2 and corresponds to the 

29 30 

Lockheed-Georgia "Wing A" ’ which was developed as part of a 
comprehensive program to acquire steady transonic experimental data 
specifically for evaluation of three-dimensional computational methods. Tne 
wing has an aspect ratio of 8.0, a taper ratio of 0.4, and leading edge sweep 
angle of 27°. A nonuniform \2% thick airfoil section having appreciable aft 
camber is employed, with approximately 5° of nose-down twist occurring 
between the root and tip sections. The section and planform geometry are 
regarded as representative of modern transport aircraft. 

Freestream conditions were established for the flow in air at M = 0.819, 

a ■ 1.96°, and Re = 8.08 x 10°. These values duplicate the test 
03 29 30 

conditions of Hinson and Burdges ’ as well as the numerical results of 
31 

Streett . Steady-state solutions were generated by the time Integration 
procedure previously summarized. In the case of the full viscous computation 
It was not possible to obtain an entirely converged steady solution due to 
boundary-layer separation at the trailing edge of the upper surface and in the 
aft cove region of the lower surface. A reasonable result was arrived at by 
the following procedure. Using the converged wedge alone solution as an 
initial profile, a converged viscous solution was generated at a Reynolds 
number of 8.08 x 10^. This value is 10 times that of the test condition and 
was sufficiently high to preclude boundary-layer separation. The computation 
was then restarted with the converged viscous result as the initial profile 
and run for 80 time steps at the test Reynolds number. In both calcuations 
the boundary-layer computation was updated each time step. The final solution 
Is probably resolved to plotting accuracy in the pressure distribution. Any 
further calculation rapidly produced separation. 



Results of the calculations are shown in Fig. 3 in terms of the wing surface 
pressure distributions where the data of Hinson and Burdges has been provided 
for. comparison. The spanwise mesh distribution was selected to coincide with 
experimental data measuring stations. Comparison of the full viscous result 
with the Inviscld solution indicates the improvement due to the boundary-layer 
calculation. It Is noted that the shock Is weakened and displaced forward by 
the Inclusion of viscous effects. Aft decambering due to boundary-layer 
displacement is quite noticable, particularly on the aft lower surface. The 
wedge alone solution Is seen to produce some anomalous behavior on the 
rearward portion of the upper surface. This is probably due to the post shock 
re-expansion and the high aft camber. The only spurious result from the full 
viscous solution is noted at n /n t ^p = 0.70 Immediately downstream of 
the shock, and may be caused by the merging of the wedge and boundary-layer 
slope modifications. In Ref. 31, the results of calculations from a pilot 
version of XTRAN3S were presented for the same case. In those calculations, 
the spike downstream of the shock shown in Fig. 3d was not observed. 


In Figure 4 the full viscous result at a 0-70 is compared with 

* *P 3? 

another computational solution. The numerical results of Streett consist 
of a full potential inviscid outer flow solution and a fully three-dimensional 
Integral boundary-layer and wake solution coupled via the displacement 
thickness. While the flow conditions for Streett's solution are the same as 
noted previously, the configuration consists of a half-body fuselage model 
with a mid-mounted wing and corresponds to another of the test cases of Hinson 
and Burdges. Streett's results compared extremely well with experiment, thus, 
the data for the wing-body combination has been omitted. Due to the presence 
of the body, the shock lies somewhat ahead of its location for the wing 
alone. The lower surface pressure distribution Is virtually the same in both 
cases. 


Upper surface distributions of the displacement and momentum defect thickness 
at midspan appear In Fig. 5. Influence of the shock can be seen near 
= 0.4. Downstream of the trailing edge, 0 and 6* tend to constant 
values which is consistent for wakes without pressure gradient. 


In the case of fully three-dimensional unsteady transonic flows over wings, 
little experimental data suitable for validation purposes is available. 
Because of this, a representative sample calculation was selected in order to 
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Illustrate results of the numerical method. For this example a section 
geometry corresponding to the M3B-A3 airfoil was chosen. The section, 
shown In Fig. 6, has a blunt leading edge, a thickness ratio of 8.9i, moderate 
aft camber, and has been selected as an AGARD standard for evaluating 
transonic aeroelastlc analysis methods. A wing planform Identical to the 
Lockheed-Georgl a "Wing A" was employed (AR ■ 8.0, x “ 0.4) except that the 
leading edge sweep angle was set to 30°. The MBB-A3 section was chosen to 
avoid the problems with boundary-layer separation evidenced In the Lockheed 
Wing A calculation. Even though the airfoil Is thinner. It still may be 
viewed as representative of a typical modern transport aircraft configuration. 

Steady-state solutions were generated for freestream conditions corresponding 
to M ■ 0.85, a s 1.0°, and Re » 10 . Results of the pressure 
distributions at half of the spanwise mesh stations are compared In Fig. 7. 
Viscous effects near shocks and on the aft lower surface are apparent. 
Corresponding chordwise distributions of the displacement and momentum defect 
thickness near the midspan station appear In Fig. 8. 

Using the steady-state solutions as Initial conditions, unsteady calculations 
were performed by specifying a forced rotation about the wing leading edge 
root such that a 3 1° + 1° sin t. The reduced frequency, k c , was 

chosen to be 0.2. Integration In time was allowed to proceed for two cycles 

of rotation. Figures 9 and 10 Indicate the instantaneous surface pressure 
distributions for the Invlscid, wedge alone, and fully viscous cases. The 
value of time here Is referenced to the beginning of the second cycle and 

given In radians. At the maximum angle of attack, a 3 2°, (Fig. 9) a 

large difference between the solutions, particularly near shocks Is evident. 
While the difference between the solutions Is less substantial at the minimum 
angle of attack, a ? 0° (Fig. 10), aft decambering due to viscous 
displacement Is still appreciable. Corresponding values of the displacement 
and momentum defect thickness are given In Figures 11 and 12. It Is noted 
that the large differences In 6* and e at the extremes of the oscillation 
cycle are brought about by the alternate formation and disappearance of the 
shock. In Fig. 13 the total wing unsteady lift and moment coefficient time 
histories are presented along with the angle of attack variation. Although 
there is little phase difference between the inviscid and viscous solution for 
this particular case, the difference in both mean value and magnitude of 
oscillation is significant. 
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VI. DISCUSSION AND CONCLUSIONS 


A method has been presented for computing the unsteady three-dimensional 
transonic flowfield about thin wings of arbitrary planform and section 
geometry Including the effects of Invlscid/vlscous Interaction. Dominant 
effects of the shock-boundary-layer Interaction were accounted for by a simple 
compuatlonal artifice, while viscous regions adjacent to the wing surface and 
along the trailing wake were described by a set of Integral equations 
appropriate for two-dimensional turbulent flows. Although this analysis was 
necessarily simplified, comparison with steady three-dimensional experimental 
data and a solution of more exact equations for a typical transport 
configuration Indicates that It will be adequate for many practical 
applications. In addition, an example calculation of a wing oscillating in 
rotation Indicated appreciable variation In unsteady aerodynamic force 
coefficients when viscous effects are accounted for. This behavior is typical 
of wings having advanced airfoil sections and may be a significant 
consideration when performing aeroelastic analysis. 

Several comments regarding the operation and possible Improvement of the 

modified version of XTRAN3S are appropriate here. No provision has been made 

In the present work for the existence of a laminar boundary layer. The reason 

for this Is that the lag entrainment equations do not perform well at low 

Reynolds numbers. Ideally, a simple attachment line model or stagnation point 

solution followed by a laminar calculation and suitable transition criteria 

could easily be Incorporated in order to establish initial conditions for the 

lag entrainment equations without altering the basic computational method. 

For the results considered here, the Reynolds numbers were sufficiently large 

such that for practical purposes the entire boundary layer could be considered 

turbulent and the upstream displacement effect could be Ignored. This may not 

be possible in all cases. Wake curvature effects have been neglected entirely 

In the present work. These may be accounted for by a simple alteration of the 

32 

basic wake jump condition, Eq. 10a . It is Indicated by the results of 
32 

Streett that wake curvature may have a significant impact upon the lower 
surface pressure distribution for wings having advanced airfoil sections. 
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The Lockheed-Georgla "Wing A" proved to be a severe test of the computational 
method. Separation was Indicated for the experimental configuration due to 
the combination of thickness, camber, and flow conditions. This evidences a 
need for treating cases of practical interest where boundary-layer separation 
occurs and requires a reformulation of the boundary-layer solution method In 
the "Inverse” mode. It should also be noted that the version of XTRAH3S 
employed for the calculations was designed to operate efficiently on the 
CRAY-IS computing system. The boundary-layer calculation method was 
Implemented as a simple modification of the basic code, but does not take 
advantage of the vectorizatlon capabilities of the CRAY-IS. This slows the 
viscous calculation somewhat, but It could easily be made more efficient. 

It is not recommended that the wedge alone be used on wings having advanced 
airfoil sections. As was noted for the cases presented here, use of the wedge 
alone tended to produce spurious post shock re-expansions. This effect may be 
reduced by varying the wedge parameters { r , r , r } but that was 

0 P K 

not done here. For wings having more conventional sections It Is expected 
that the wedge alone will be quite useful. 

The method of time Integration of fluid dynamic equations of motion offers a 

unique capabilty for performing aeroelastic analysis. Unsteady structural 

equations of motion governing the deformation of an elastic wing may be 

coupled with these directly and simultaneously Integrated in time to determine 

•the aeroelastic response of the combined system. This technique has already 

13 

been applied for a purely invlscid flowfield, and may easily be extended 
to include viscous effects by the described method. Addition of the Integral 
equations for viscous calculations increases computation time (CPU) by about 
20S, thus maintaining the basic efficiency of the'invlscld algorithm. This 
consideration provided much of the motivation for employing a simplified 
viscous analysis, and appears to be justified for the purpose of structural 
applications. 


31 



APPENDIX A 


Conditions at the edge of the boundary layer may be obtained from the invlscld 
solution along the wing surface or wake centerline evaluated on 5 a 0*. 

In particular, using the small disturbance assumption and simple sweep theory, 
we may write 


H e »{l * [1 ♦ M, 2 ] <i> jn, 

p e _ . /It u 2. 

5T * 1 - 'c 1 V 


COS A s , 


- 1 - ( Y - 1) M^ 2 (I) 4» r 


( A-l ) 
(A-2) 

(A-3) 


The viscosity ratio Is obtained from the Sutherland viscosity law In the form 


CO 


e%3/2 


S o 

1 +T“ 

o 

7~ 

L’p + ’pJ 

CO CO 


(A-4) 


where S Q Is a constant for the specific gas under consideration 

(S Q = 110°K. for air). The adiabatic recovery factor, r, is chosen as a 

function of the turbulent Prandtl number 


r 



(A-5) 


The following parameters then completely define the lag entrainment equations 


1 o 1/2 

F c ■ Cl + P-J-i) r fl fi 2 ], 

F r a 1 + 0.056 M e 2 , 


Re 0 


p e U e 9 

^e 


(A-6) 

(A-7) 

(A-8) 
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C fo“ 


<r» “ + ^ cos A s £ ^ £• 

0.01013 _ 0.00075] for wall boundary layers 

(A-9) 
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C 
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C- 


0.9 


(— - 0.4) 


0.5] C fo for wall boundary layers 


Ho 


l 0 for wake shear layers, 

H q - (H + 1) [1 + (^-) r M fi 2 ] - 1 


dH 


- - [ (*-!)* ] 
^1 1.72 + 0.02 (H - 1)^ 
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1 H- 1 


F A° 


0.02C^ + C E 2 + ( 0, ° C fo ) 
3 
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C t * (1 + 0.1 M e 2 ) (0.024 C £ + 1.2 C E + 0.32 C fQ ), 
rl for wall boundary layers 


Ll/2 for wake shear layers, 
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dU 


e c- f 


\ ^ EO- "0 
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APPENDIX B 

MODIFICATION OF XTRAN3S INPUT DATA REQUIREMENTS 
FOR YISCOUS CALCULATIONS 


Tho input data required p, XTRAN3S Is In the form of a card deck or a card 
trnaje file. All data Is free field. Section A describes the Input data deck 
tructure while Section a spells out the reared Input cards for 
mplementation of the viscous option and the associated formats and 
groundrules. This Appendix may be considered as a revision of Section V of 
Ref. 26, and should be used in conjunction with that Section. 

A. Input Data Deck Structure 

The input data deck for XTRAN3S has four levels of organization. 

1. Program Deck - This includes all Inputs (other than control cards which 
are discussed ,n Section VIII of Ref. 26) for one problem. The deck 
boundaries are 


first record - BEGIN PROBLEM DEFINITION 
last record - END OF PROBLEM 

2. Data Section - A program deck is divided Into data sections. 

Specifically, the deck Is divided Into the following ten data sections 

1) Problem Definition Section 

2) Computational Control Section 

3) Computational Grid Section 

4) Geometry Section 

5) Boundary Condition Section 

6) Structural Modal Section 

7) Structural Matrix Section 

8) Checkpoint/Restart Section 
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9) Post-processing Section 

10) Viscous Calculation Section 

As noted previously, the first section must be the ProDlem Definition 
Section. The order of the remaining sections Is immaterial. The first 
record In each section must be of the form "BEGIN etc." No section 
terminator Is required. 

Data Record - Sections are In turn divided Into data records or 
statements. As noted above, the first record is of the form "BEGIN 
etc." In most cases the order of the remaining records in the section 
Is Immaterial. All exceptions to this rule are specifically noted. 
Record boundaries are governed by the following rules. 

o If the last non-blank character on the card Isa "+" then the 
record continues onto the next card. 

o Maximum record length is 250 items. 

o Record terminators are / or card boundaries. 

o The space between / and the card boundary is "ignored" and hence is 
available for comments. If a record is only comment it must begin 
with */. 

Data Item - Records are in turn composed of data items. This Is the 
finest level of subdivision for the program deck. The delimiters or 
"boundaries" of items are of two types. 

o colon : this is used at most once per record and then only after 
keywords. 

o commas, blank spaces and record terminators. 


■? myu y ujj ,-upjh 







4! . 


^ijnii^ii.i.. mjp ii um»ii il 
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B. Record Formats 

The following ground rules concerning notation apply throughout this section. 

a) Ct ]]— Data items enclosed by double brackets have optional input 
formats. One or more of the Indicated options must be selected. 

b) [ ]-- Data items enclosed by brackets have optional Input formats. 

Only one of the Indicated options must be selected. 




L ! 
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c) ( )— Data items enclosed by parentheses have default values. If the 
default Is acceptable for definition of the problem data, the particular 
Item or items need not be Input. All default values are defined in the 
descriptions of the input data. 

d) .ITEM — An item typed In all upper case letters is called a key-word. 

At least the underlined portion of a key-word must be input. This Is 
always the first four characters Including trailing blanks. 

e) Item — An item with only the leading character typed in upper case 
denotes that it must either be selected from a list of system key-words 
or that it is Identical to an Item previously defined by the user. 

f) item — An item typed In all lower case letters is defined strictly by 
the user. 

The formats for data sections 1 through 9 are given in Ref. 26. The 
additional data required for viscous calculations Is given in Section 10, 
which follows. 

10. Viscous Calculation Section 

The purpose of this section Is to specify the Viscous Calculation option 
if the user so desires. If this option is not desired, then this 
section should be omitted. 
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1) BEGIN VISCO US FLOW PARAMETERS 

This data record is used to indicate that processing of data 
associated with the Viscous Calculation Is to follow. The 
following card statements (2-5) are optional and need only be input 
If values other than the defaults are to be specified. Data 
records 2-4 define data associated with the viscous wedge (all 
values are nondimensional by local section chord). 

2) ( SHOCK OFFSET DISTANCE: xoffst) 

Input of the shock offset distance. 

The default is SHOCK OFFSET DISTANCE: .02. 

3) (PRECURSOR LENGTH: xprec) 

Input of the precursor length. 

The default is PRECURSOR LENGTH: .02. 

4) ( RAMP LENGTH: xramp) 

Input of the ramp length. 

The default is RAMP LENGTH: 0.10. 

5) (CALCULATION INTERVAL: iblcal) 


Input of the viscous fiow calculation update Interval for both 
wedge and boundary-layer computations. 

The default is CALCULATION INTERVAL: 1 

6) BOUNDARY-LAYER SOLUTION 

This data record is required if the boundary-layer calculation Is 
to be executed and it also is the first data record associated with 
the boundary-layer calculation. The following card statements 
(7-11) are optional and need only be input if values other than the 
default values are to be input. 


. ". . X .. ■ ' •' ^ ' . j 


7) ( PRINT INTERVAL: Iblprt) 

Input of the boundary-layer print interval. 
The default is PRINT INTERVAL: 50. 

8) ( REYNO LDS NUMBER: reyinf) 


Input of the free stream Reynolds number based on root chord. 
The default is REYNOLDS NUMBER: 1.0E7 


9) ( TEMPE RATURE: tinf) 


Input of che free stream temperature (degrees Kelvin). 

The default is TEMPERATURE: 300.0 

10) ( SUTHE RLAND LAW CONSTANT: SO) 

Input of the Sutherland Law Constant (degrees Kelvin). 

The default (defined for air) is SUTHERLAND LAW CONSTANT: 110.0. 

11) ( PRAND TL NUMBER: prt) 

Input of the turbulent Prandtl number. 

The default is PRANDTL NUMBER: 0.9 


I 

i 

I 


1 
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■ APPENDIX C 


The n-roesh distribution and wing surface slope definition for the two 
geometries considered are given here. In each case the wing surface is 
represented as an analytic function having the following form: 

+ a i/2? 1/2 + a ! S 1 + a 2 E 2 + a 3 E 3 + a^ 4 > a^ 5 . 

1. Lockheed Georgia “Wing A" 


a. n -mesh 


Index 


1 

- 0.21 

2 

0.00 

3 

0.21 

4 

0.42 

5 

0.63 

6 

0.83 

7 

1.11 

8 

1.39 

9 

1.94 

10 

2.25 

11 

2.50 

12 

2.65 

13 

2.75 

14 

2.80 

15 

2.90 

16 

3.05 

17 

3.30 

18 

3.80 

19 

4.30 

20 

5.30 
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b. 

surface geometry 





Root 


li& 



Upper Surface 

Lower Surface 

Upper Surface 

Lower Surface 

a o 

0.023996815 

0.023517482 

-0.017931003 

-0.017921293 

a l/2 

0.210881410 

-0.148954440 

0.175083370 

-0.192270330 

a l 

-0.348816530 

-0.010928663 

0.019118126 

0.371830370 

«, 

0.418018550 

-0.106489380 

-0.494377680 

-0.938122680 

a 3 

-0.637010050 

0.391707720 

1.12296990 

1.757023200 

a 4 

0.395787680 

0.150884270 

-1.211785000 

-0.989446110 


-0.085424090 

-0.324138520 

0.426635820 

0.024627938 
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2 
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6 
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10 

2.25 

11 

2.50 

12 

2.65 

13 

2.75 

14 

2.80 

15 

2.90 

16 

3.05 

17 

3.30 

18 

3.80 

19 

4.30 

20 

5.30 
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Surface Geometry 


Root 


HE 


Upper Surface 

Lower Surface 

Upper Surface 

Lower Surface 

0.0 

0.0 

0.0 

0.0 

0.1064264 

-0.1064262 

0.1064264 

-0.1064262 

0.0051279 

0.2190702 

0.0051279 

0.2190702 

-0.1236988 

-0.6049436 

-0.1236988 

-0.6049436 

-0.0149070 

0.9704346 

0.0149070 

0.9704346 

0.0270518 

-0.4781350 

0.0270518 

-0.4781350 

0.0 

0.0 

0.0 

0.0 
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Figure 1. Viscous Ramp Geometry 
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Figum 6. MBB-A3 Airfoil Section 
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Figure 7a. Steady Surface Pressure Distributions for MBB-A3 Wing at G/Hf/p = 0.00 
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Figure 7 b. Steady Surface Pressure Distributions for MBB-A3 Wing at T)/T) 
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Figure 7 c. Steady Surface Pressure Distributions for MBB-A3 Wing at F|/r )tip ~ 0.36 
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Figure 7d. Steady Surface Pressure Distributions for MB 8- A3 Wing at T) /T) tlf} =* C.6J 
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